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The amplitude and phase of Shubnikov-de Haas oscillations have been analyzed in detail for 
two-dimensional electron gases subjected to a weak unidirectional periodic potential modula- 
tion. The amplitude is suppressed, accompanied by inversion of the phase, at the maximum 
bandwidth conditions at low magnetic fields. The suppression is gradually taken over by the 
enhancement with the increase of the magnetic fields. The suppression and the enhancement 
are attributable to the coUisional and the diffusion contribution of the modulated potential to 
the conductivity, respectively, the former (the latter) being dominant at low (high) magnetic 
fields. A theoretical calculation that takes the two types of contributions into account shows 
semi-quantitative agreement with experimental traces. 
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1. Introduction 

The Shubnikov-de Haas oscillation (SdHO) in a two- 
dimensional electron gas (2DEG) is a manifestation of 
Landau quantization, En = {N -\- l/2)hujc, in the mag- 
netoresistance."'^ The oscillation is periodic in and 
damps with decreasing magnetic fields since the cy- 
clotron energy hwc = heB/m* (with m* the electron 
effective mass) diminishes with respect to both the ther- 
mal blurring of the Fermi energy and the dis- 
order broadening F of the Landau levels (LLs). In a 
unidirectional lateral superlattice (ULSL), a 2DEG sub- 
jected to a one-dimensional (ID) periodic modulation, 
another analogous oscillation periodic in 1/B emerges. 
The oscillation, known as the commensurability oscilla- 
tion (CO),^"** results from the commensurability between 
the cyclotron radius Rc = hk-p/eB and the period a of 
the modulation, and also damps with decreasing B. Here 
fcp = represents the Fermi wave number with 

the areal density of electrons. In the CO, the energy dif- 
ference between adjacent fiat-band conditions (see be- 
low), {akF/2)hujc, takes the place of hu!c in the SdHO. 
Since {ak-p/2) 3> 1 in most experiments performed so 
far, the SdHO vanishes and only the CO survives at the 
low end of magnetic fields. At low enough temperatures 
where thermal damping is not so severe, however, the two 
classes of oscillation coexist over a range of magnetic field 
with an intriguing interplay. The behavior of the SdHO 
under the infiuence of the simultaneously present CO is 
the main subject of the present study. 

A unidirectional periodic modulation of the electro- 
static potential. 



V{x) — Vq cos{qx) 



(1) 



with q = 27r/a, lifts the degeneracy of the LLs. The en- 
ergy becomes dependent on the position xq of the guiding 
center and reads, in the first order perturbation theory 
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valid for a weak modulation Vq <ti Ep, 

EN,qxo = (^^ + ^^'^ + ^^'B cos{qxo), (2) 

where Vn,b = Voe^'^/^Li^{u) with Ln{u) the Laguerre 
polynomial, u — q'^P/2, and / — y/h/eB the magnetic 
length. Vn,b at the Fermi energy can be approximated, 
using an asymptotic expression for A'^ ^ 1 appropri- 
ate for low-magnetic-field range relevant to the present 
study, by. 



Vb^Vo 



nqRc 



■ cos 



{qRc 



(3) 



Vb oscillates with B, and the width of the Landau bands 
2|yB| takes maximum at 

2i?c , 1 
a 4 
and vanishes at the flat band conditions 



in= 1,2,5,...), 



(4) 



2Rr. 



1 

= n 

4 



(n= 1,2,3,...). 



(5) 



The oscillation of the Landau bandwidth is the origin 
of the CO and, at the same time, is responsible for the 
modulation of the amplitude and the phase of the SdHO. 

The modulation of the SdHO amplitude was first re- 
ported by Overend et al.^ They exploited a periodic mod- 
ulation of magnetic field instead of electrostatic poten- 
tial. For a magnetic field modulation, eqs. (4) and (5) 
interchange their roles, with eq. (4) representing the fiat 
band conditions. The underlying physics, however, is ba- 
sically the same for both types of modulations. The au- 
thors reported that the SdHO amplitude remains large 
at the fiat band conditions [eq. (4)], while is suppressed 
at the maximum bandwidth conditions [eq. (5)]. Fur- 
ther study by the same group^' ^ showed the phase inver- 
sion of the SdHO at the maximum bandwidth conditions 
when the bandwidth is larger than hujc. They explained 
their observation as an effect of the modulated density 
of states (DOS), which affects the conductivity through 



2 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



the collisional contribution. The phase inversion is essen- 
tiaUy the same phenomenon as the even-odd fiUing-factor 
switching reported for a ULSL with strong electrostatic 
potential modulation in the quantum Hall regime,* and 
is attributed to the van Hove singularities in the DOS. 

Similar suppression of the SdHO amplitude at the 
maximum bandwidth conditions [eq. (4)] was also re- 
ported for electrostatic ULSLs with a strong modula- 
tion amplitude induced by a patterned InGaAs stressor 
layer. ^ Interestingly, the authors observed that the trend 
is reversed at higher magnetic fields; the amplitude is 
enhanced at the maximum bandwidth conditions. The 
origin of this enhancement, as well as of the inversion of 
the tendency with the magnetic field, still lacks compre- 
hensive explanation. 

The purpose of the present paper is to achieve more 
quantitative understanding of the behavior of the SdHO 
under periodic modulation. In the previous works quoted 
above, ULSLs with a large amplitude periodic modu- 
lation are employed, which is obviously advantageous in 
attaining the modulation of the SdHO strong enough to 
be readily observed. However, the amplitude of the peri- 
odic modulation or the width of Landau bands exceeding 
10% of the Fermi energy seems to be rather incompati- 
ble with the perturbative treatment as in eq. (2). In the 
present work, we keep IVbI/E'f < Vq/Ef to be less than 
five percent, which validates the comparison of the ex- 
perimental data with perturbation theories. We employ 
low temperatures, which allows us to observe the SdHO 
down to low magnetic fields (~ 0.05 T) where Vo/fiuJc 
becomes large even for our small Vq. 

2. Experimental 

We examine two ULSL samples with slightly differ- 
ent characteristics, as tabulated in Table I. The two 
samples are fabricated from two GaAs/AlGaAs 2DEG 
wafers having nominally the same structure^" but dif- 
fering in the carrier mobilities owing to the conditions of 
the molecular beam epitaxy (MBE) chamber used for the 
growth. The electrostatic potential modulation is intro- 
duced via strain-induced piezoelectric effect by placing 
negative electron-beam (EB) resist on the surface. As 
depicted in the right inset of Fig. 1 (a), each sample is 
patterned into a Hall bar that has two sets of voltage 
probes to measure the longitudinal and Hall resistivity 
of the section with (ULSL) and without (plain 2DEG) 
the potential modulation at the same time. The values 
of the mobility /i and the density tabulated in Table 
I are those measured after brief illumination by LED, 
the condition in which the magnetoresistance traces pre- 
sented in this paper are taken. No sign of deterioration 
in fi or change in Ue arising from the microfabrication 
process is discerned. 

For Sample B, the periodic modulation is introduced 
by a EB-resist grating having a conventional periodic 
line-and-space pattern. In this approach, higher harmon- 
ics inevitably mixes in the potential profile as the pe- 
riod a becomes large compared to the depth d {— 90 
nm for the present samples) of the 2DEG plane from 
the surface. An unconventional strategy of employing 
a quasiperiodic pattern is applied for Sample A with 



Table I. Sample parameters 





Sample A 


Sample B 


Period a (nm) 


231 


184 


Amplitude Vq (meV) 


0.39 


0.31 


Mobility fj, (m^/Vs) 


69 


101 


Quantum mobility fiQ (m^/Vs) 


11.9 


7.2 


Density rie {10^^ m~^) 


2.9 


2.9 




S(T) 



Fig. 1. (a) (Color online) Magnetoresistance traces for the ULSL 
(solid line) and the plain 2DEG (dotted line) taken at 15 mK 
for sample A. A trace for the ULSL at 4.2 K is also shown with 
dot-dashed line.^^ The upper right inset depicts the schematic 
configuration of the sample, (b) The absolute value of Vb [eq. 
(3)], with Vb > (< 0) plotted by solid (dashed) line. 



a larger period. Here, slabs of EB-resist are placed on 
the "L"s of the Fibonacci sequence "LSLLSLSL..J' , 
with L=104 nm and S'=64 nm [L/S is set to 0=(1 + 
\/5)/2=1.618...]. Generally in such Fibonacci ULSLs, the 
analysis of the CO reveals several frequency components, 
each corresponding to one of the self-similar generations 
of a potential profile mutually scaled by the factor cj)}^' 
In the particular sample explored in this study (Sample 
A), a single component with effective period a=231 nm, 
corresponding to the average distance between adjacent 
"S"'s, overwhelmingly dominates the potential profile, al- 
lowing the profile to be virtually regarded as periodic. 
Although somewhat counterintuitive, a simple sinusoidal 
potential profile described by eq. (1) is realized better in 
Sample A than in Sample B because of the absence of 
the higher harmonics. 

Cursory characterization of the ULSL samples can 
be done by performing a Fourier transform of the 
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Fig. 2. (Color online) Fourier transform of (d? /dB'^)[pxx{B)/ po] 
versus l/B for the ULSL in sample A (a) and B (b). Fourier 
spectra for the plain 2DEG is also shown for (a) by dotted line 
with the shade underneath. Peak positions are given by combi- 
nations of integer multiples of /s^h and fco, where fgdH (/co) 
represents the fundamental frequency for the SdHO (CO). 

(d'^ / dB^)[pxx{B) / pq] vs. 1/B curve, where pxx{B) rep- 
resents longitudinal magnetoresistance and po = Pxx{Q)- 
The second derivative conveniently eliminates the slowly- 
varying background from the magnetoresistance traces. 
The Fourier spectra shown in Fig. 2 exhibit peaks corre- 
sponding to the CO (/co), the SdHO (/sdn), and their 
combinations (p/sdni'z/co with p, q integers). The CO 
peaks up to the fourth harmonic (/co, 2/co, 3/co, and 
4/co) are seen for Sample B, indicating the presence 
of the corresponding harmonic contents in the potential 
profile. For Sample A, in contrast, only one CO peak is 
found, justifying the description by eq. (1) of the poten- 
tial profile. The presence of the peaks at p/sdH ± <lf GO 
reveals the interplay between the two types of oscilla- 
tion that will be discussed in detail in the following sec- 
tions. The higher harmonics of the CO in Sample B do 
affect the interplay as indicated by the presence of a peak 

/sdH — 2/co-^^ 

More quantitative account of the actual potential pro- 
file requires analyses of the CO amplitude. The values of 
the fundamental component Vb tabulated in Table I are 
obtained by such analyses that take the effect of damping 
by scattering into account. Higher harmonic contents 



can also be evaluated by analyses using Fourier band pass 
filters, and the second, the third, and the fourth har- 
monic components for Sample B are found to be V2=0.10 
meV, ^3=0. 07 meV, and ^4=0.05 meV, respectively. Al- 
though these higher harmonics can in principle compli- 
cate the analysis below for Sample B, we neglect them 
for simplicity and assume, in the rest of this paper, that 
both Sample A and Sample B have a potential profile of 
the form given by eq. (1). 

Measurements are performed in a top-loading dilution 
fridge at the base temperature T ~ 15 mK. Standard 
low-frequency (13 Hz) ac lock- in technique with a cur- 
rent /ims^lO nA is employed for the resistivity measure- 
ment. We have checked that the electron heating by the 
current is negligible in the magnetic-field range of the 
present interest {\B\ < 1 T) by reducing /ims down to 
0.5 nA, for which the magnetoresistance trace shows no 
difference except for much worse signal-to-noise ratio. A 
low sweep rate dB /dt=O.Ol T/min is employed in or- 
der to avoid undesired hysteresis of the superconducting 
magnet. The remnant hysteresis is further calibrated by 
using the simultaneously measured Hall resistivity. 

3. Results and Discussion 

3.1 Experimentally obtained Shubnikov-de Haas oscil- 
lation 

Figure 1 (a) shows magnetoresistance traces of Sam- 
ple A for both the ULSL and the adjacent plain 2DEG. 
Magnetoresistance of the ULSL at 4.2 K is also shown, ^"^ 
which essentially represents the pure CO for \B\ < ^0.5 
T where the SdHO has already damped out. In the low 
magnetic field range {\B\ < ~ 0.25 T), the SdHO of the 
ULSL is suppressed at the maxima of the CO, while re- 
mains almost unaltered from that of the plain 2DEG at 
the minima of the CO. In contrast, the SdHO amplitude 
is observed to be enhanced at the CO maxima for higher 
magnetic fields. This is most clearly seen for the CO peak 
at 0.64 T. In Fig. 1 (b), the half width of Landau bands 
calculated by eq. (3) are plotted. It can readily be con- 
firmed that, as is well known, "^''^ the maxima and the 
minima of the CO correspond to the maximum band- 
width and the fiat band conditions, respectively. 

To look into the details of the behavior of the SdHO, 
the rapidly oscillating parts of the magnetoresistance, 
^Psdn/po, are extracted and plotted against the Landau 
level filling factor v=ne,h/ eB in Fig. 3 (a) for both the 
ULSL and the plain 2DEG of Sample A. The correspond- 
ing magnetic field is shown on the top axis. The extrac- 
tion of the Apsdn/po is done by applying a Fourier high 
pass filter to the Pxx/ Po vs. 1/B curve, with the threshold 
set at a frequency higher than /co, and further by sub- 
tracting the average of the upper and lower envelopes, 
as was done in the analysis of the CO."'^^ The SdHO 
is known to damp with decreasing magnetic field by 
the thermal damping factor A{T/Tc) and the scattering 
damping factor exp(— tt/wcTq), where A{x)=x/ sinh(x), 
fcBT'c=^c/27r^, and tq the single particle (quantum) 
scattering time. Traces shown in Fig. 3 are Apsdn/po di- 
vided by these damping factors, applying to both traces 
the value of the quantum mobility pq — eTQ/m* — \l.% 
m^/Vs determined from the actual damping of the SdHO 
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Fig. 3. (Color online) (a) Experimentally obtained SdHO, divided by the thermal and the scattering damping factors, for sample A. 
See text for details, (b) Density of states divided by the exponential damping factor, calculated by eq. (7) or eq. (10) with Vq = 0.49 
meV. The thin line with the shade underneath and the thick line represent the plain 2DEG and the ULSL, respectively, plotted against 
the filling factors (bottom axes) or magnetic fields (top axes). The half width of the Landau bands, \Vb\, is also plotted in (b) (right 
axis). 



of the plain 2DEG. The factor A{T/Tc) turns out to 
barely deviate from unity, reflecting the fact that the 
thermal damping is negligibly small at this low temper- 
ature. Since spins are unresolved for the magnetic-field 
range in the present study, the minima and maxima of 
the SdHO are expected to take place at even and odd fill- 
ing factors, respectively. This is actually what we observe 
for the plain 2DEG. For the ULSL, the SdHO amplitudes 
are suppressed at the maximum bandwidth conditions, 
B=0.239, 0.183, 0.148, 0.125, 0.107, 0.094, 0.084, 0.076, 
0.069, 0.064 T [n=3, 4,.. .,12 in eq. (4), note the \Vb\ plot- 
ted in Fig. 3 (b)] . The SdHO amplitude at the suppressed 
conditions decrease with decreasing magnetic fields until 
it vanishes at B ~ 0.125 T, and then revives at still lower 
magnetic fields but with the position of peaks and dips 
inverted; there, the minima and maxima occurs at odd 
and even filling factors, respectively. This is the even-odd 
transition reported in refs. 6-9, which was attributed, by 
using numerically calculated DOS and the resultant con- 
ductivity, to the broadening and the van Hove singularity 
of individual LLs that result in the peaks of DOS at even 
filling factors. 

3. 2 Comparison with calculated density of states 

It is well established that the ApsdH/[poA{T/Tc)] is 
proportional to the oscillatory part of the DOS at the 
Fermi energy, AD{Ep)/Dq, at low temperatures pro- 
vided the |ApsdH/Po^l is not too large. ^'^^ In this sub- 
section, we make a detailed comparison of the experimen- 
tally obtained SdHO described in the previous subsection 
with a calculated DOS. 



First we deduce an analytic expression for the DOS 
under a periodic modulation eq. (1), instead of resort- 
ing to a numerical calculation as was done in the pre- 
vious works. We start by recalling the DOS for a plain 
2DEG.^ The disorder broadened line shape of each LL 
peak is approximated here for simplicity by a Lorentzian, 
P{E) = (r/7r)/(^2 _^r2), with the width F independent 
of B. The DOS is obtained by summing up the LL peaks, 
including the factor 2 for spin degeneracy, 



D{E) 



Y,P{E-Em) 



N=0 



2E 



cos 



fc=i 



27rfc e 



(6) 

where Dq = m* /nh? represents the constant DOS of a 
2DEG in the absence of magnetic field. We introduced 
dimensionless parameters, e = E/fiiOc and 7 = T/fiWc- 
In the second line of eq. (6), we made use of the Poisson 
sum formula. Since exp(— 27r7) <gC 1 for small magnetic 
fields, it is usually a good approximation to keep only 
the k=l term in the summation: D{E) ~ Dq + ADi{E) 
with 

ADi{E) 



Do 



= —2 cos(27re) cxp(— 27r7). 



(7) 



The proportionality Apsdn/po^ oc ADi/Do imphes 
F = H/2tq by the comparison of the exponential fac- 
tors. Therefore the quantum mobility fiQ = 11.9 m^/Vs 
corresponds to F = 0.073 meV. 
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Fig. 4. Plots of Jo{2'Kx) (a) and Ji{2'kx)/{2-kx) (b). 



In Fig. 3 (b), we plot the oscillatory part of the 

DOS at Ey = fic/Do^'^ divided by the damping 
factor, ADi(£;F)/[£'oexp(-27r7)] = -2cos(cf) with 
£f = E-p/tkjJc = i'/2. Comparison of Figs. 3 (a) and 
(b) confirms that the relation ApsdH/[/Oo^(T/Tc)] cx 
ADi{Ef) / Do actually holds for our plain 2DEG, aside 
from the deviation at higher magnetic-field range |-B| > 
0.4 T where the peak height of ApsdH/[paA{T/Tc)] di- 
minishes accompanying the onset of spin-splitting, which 
was not taken into consideration in the calculation of the 
DOS. 

Upon switching on the potential modulation eq. (1), 
each LL peak further broadens from P{E) to (with t = 
QXo) 



P{E, Vb) = - [ P{E - Vb cost)dt. 

Jo 

This can be formally rewritten as 
±( I 

2w \y/E-VB+ iVy/E + VB+iT 



(8) 



y/E-Vn- iTVE + Vb- iT , 

and displays broadening or, when Vb is large compared 
to F, splitting (van Hove singularities) of a LL peak. 
Correspondingly, the DOS becomes 



D{E,Vb) 



N=0 



P{E-En,Vb) 



^ 1 r 
; > — / cos 



2nk [e — vb cos t~ 2 



dt 



xe 



-27rfc7 



2 ^ cos 
fe=i 



27rfc 



Jo {2-KkVB) 6-2'^'=^ 



(9) 

with Vb = Vs/fioJc, and Jo{x) the Bessel function of 
order zero. ADi{E) acquires an extra factor Jo(27rt;B): 

ADi (E) 

= -2Jo{2nVB) cos(27r£) exp(-27r7). (10) 

With the decrease of B, vb oscillates periodically in 
1/B back and forth around vb — 0, increasing its ampli- 
tude proportionally to B~^^^ [see eq. (3)]. As shown in 
Fig. 4 (a), the function Jo(2nVB) decrease from 1 with 
the increase of \vb\, becomes zero at \vb\ = 0.3827... 



~ 3/8 and then changes its sign. Therefore the oscilla- 
tion of ADi takes minimum amplitude at the maximum 
bandwidth conditions while \vb \ stays less than 3/8, dis- 
appears when a maximum of \vb\ touches ^3/8, reap- 
pears with inverted sign for \vb\ larger than 3/8.^^ Thus 
the position where the oscillation of ApsdH/[poA{T/Tc)] 



vanishes is the landmark of \Vb\ 



0.3827fia)c, on 



an assumption that the relation Apsdii/[poA{T/Tc)] oc 
ADi{Ef) / Do holds also for ULSLs. From the position 
wc can determine the modulation amplitude Vq using 
eq. (3). The oscillation disappears and a small sign of 
inverted peak is observed a.t B = 0.125 T in Fig. 3 (a). 
From this wc obtain Vq = 0.49 meV. The DOS calculated 
using Vb = 0.49 meV plotted in Fig. 3 (b) reproduce the 
line shape of the observed Apsdn/ipo^iT /Tc)] in Fig. 3 
(a) quite well for B <~ 0.25 T, confirming the propor- 
tionality relation. The value Vq = 0.49 meV is roughly 
25% larger than that deduced from the CO, the reason 
of this discrepancy remains to be elucidated. 

For Sample B, the modulation amplitude is not large 
enough to achieve the collapse or the peak/dip inversion 
of the SdHO. However, the modulated SdHO trace simi- 
lar to Fig. 3 (a) is well reproduced for B <~ 0.2 T by eq. 
(10) with Vb = 0.38 meV, again some 20% larger than 
the value determined from the CO. 

For higher magnetic fields, SdHO amplitude is en- 
hanced at the maximum bandwidth conditions, as men- 
tioned earlier. This is obviously not reproduced in the 
calculated DOS and requires an alternative explanation. 
We note in passing that the onset of the spin-splitting 
for the ULSL shifts to a higher magnetic field due to the 
detrimental effect of the modulation-induced broadening 
of LLs on the Zeeman gap, and is outside the magnetic- 
field range for the present study. 

3.3 Comparison with calculated conductivities 
3.3.1 Analytic expressions for conductivities at low 
temperatures 

It has been pointed out by Peeters and Vasilopoulos 
that a periodic potential modulation V{x) alters the con- 
ductivity of a 2DEG via two different routes. The colli- 
sional (hopping) contribution corresponds to the effect of 
the DOS we have just discussed above. Peaks in the DOS 
boost the conductivity hence the resistivity, through the 
increase in the scattering rate. The other route, the diffu- 
sion (hand) contribution, results from the drift velocity, 

_ 1 ^£^Ar.(-g;^fcy) qVB 
ft 



dk,, 



eB 



sin((jf.Xo), (11) 



which enhances the ayy hence the pxx- Here, we made use 
of the relation xo = —Pky. It is this effect that mainly 
contributes to the CO, with Vy being maximum (zero) 
at the maximum-bandwidth (flat-band) conditions. It is 
worth pointing out that the diffusion contribution has no 
counterpart in a plain 2DEG and therefore vanishes with 
Vq ^ 0, while the coUisional contribution is basically an 
ordinary SdH effect with a due modification introduced 
by the periodic modulation. 

Shi et al. calculated by a perturbation theory the coUi- 
sional (c7™' = (7y°') and the diffusion conductivities (f ) 
for 2DEGs under ID modulation of both an electrostatic 
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potential and a magnetic fieldJ Taking only the elec- 
trostatic potential, eq. (1), relevant to the present study 
into account, they read. 



f 

Jo 



dtP\E-EN,t), 



(12) 



and 



^dif 



V2 
•^0 



hL0c{akF/2)' 



:2tt 



dE --4- 



dl 
dE 



°° r t2 /■'^ 

X V e-"/2Ljv(w) / dtP{E - EN,t)sm'^t, 



(13) 



respectively, where f{E) = {1 + exp[(i; - £;f)/A;bT]}-i 
is the Fermi-Dirac distribution function and ctq = 
riee^T/m* = (e'^ /h){Ep/To) represents the conductivity 
at zero magnetic field with r = h/2TQ the momentum 
relaxation time.^^ In eq. (12), contributions from inter- 
Landau-band hoppings arc omitted. The authors of ref. 
7 made comparison with experimental results by numer- 
ically evaluating these equations. Here, we will take a 
further step and deduce from cqs. (12) and (13) approxi- 
mate analytic formulae appropriate for comparison with 
our experimental data. 

Firstly, since we employ a low temperature (T ~ 15 
mK) in our measurement, the derivative of the Fermi- 
Dirac distribution function can safely be approximated 
by the delta function 6{E — Ep), resulting in, 



^col 



J2 (2iV + 1) / dtP\Ep - EM,t) (14) 



and 



^dif 



N=0 



[e--/^LM{u)f 



[ dtP{E] 
Jo 



F - EN,t) sin^ t. 



(15) 



Equation (14) can be rewritten along the same line as in 
eq. (6). Leaving the detail of the derivation to the Ap- 
pendix, the calculation leads to, up to the leading term, 



r, 



2 

7 



1 + 2 ^ (1 + 27rA;7) 



k=l 



X COS 



27rfc I ep 



Jo(27rtoB)e-2'^'=T 



(16) 

In eq. (15) the term e~"/^Z/jv(w) may be replaced by its 
asymptotic expression at the Fermi energy, since P{Ef — 
Em^i) takes significant value only at ~ A^f = [E-p/^c\ 
(the integer part of E-p/hoJc), the index of LL in which 
the Fermi level resides (A^f ^ 1 for low magnetic fields). 



and therefore 



^dif 



2TT^Vi 



huj^{akT/2f 



cos 



X - / dt sin^ t y P{Ef -ENt)- (17) 

As wc have done in eq. (9), we employ eq. (6) to rewrite 
the summation Y1^=qP{.^'p ~ Ejs!,t)^ and use relations 
{1 / tt) sm^ tdt = 1/2, {1 / tt) cos{x cost) siv? tdt = 
Ji{x)/x, and (I/tt) s,m{x cost) siv? tdt = 0, to finally 
obtain 



^dif 

yy 



T/2 
•^0 



fkVcEF{akF/2) 



cos 



(*-i) 



1+ 



4E 



cos 



k=l 



2-Kk ef 

' 2 



Jl(27rfc^B) ^-2^fc^ 
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where Ji{x) is the Bessel function of order one. The 
conductivities are translated to resistivities by the inver- 
sion of the conductivity tensor. For not too small mag- 
netic fields [B >-0.05 T), p™VPo - (wct)V™Vcto and 
Pix/Po — {'^cT)'^(^yy /cro to a good approximation. The 
resultant resistivities are 



rol 



P_ 

Po 




l-|-2^(l-F27rA;7) 



k=l 



X COS 



27rfc £f 



Jo(27rtoB)e-2'''=T 



(19) 



and 
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Jl(2TTkvB)^_27.kj 
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(20) 

The B dependence of eq. (19) inherits that of eq. (9), as 
expected, with minor discrepancy, to be discussed in the 
following subsection, resulting from the factor (l + 2Trkj) 
in the summation. The first term in eq. (20) describes the 
CO and is exactly the same as the expression given in 
ref. 23. The amplitude of the CO increases linearly with 
B, as revealed by the factor lvcT. (The damping factor 
due to scattering^^ that results in additional B depen- 
dence is not included here.) The second term represents 
the diffusion contribution to the SdHO. The amplitude 
of the oscillation is modulated by the CO, and is there- 
fore enhanced (suppressed) at the maximum bandwidth 
(fiat band) conditions; the phase of the modulation is at 
odds with that of the collisional contribution. Owing to 
the _B-linear dependence of the CO mentioned above, the 
diffusion contribution to the SdHO raise its relative im- 
portance with B and can outweigh the collisional contri- 
bution above a certain B. This qualitatively explains the 
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Fig. 5. (Color online) Experimentally obtained SdHO versus 1/B 
for sample A (left axis) and the half width of Landau bands (right 
axis). 

observed transition, with the increase of B, from the sup- 
pression to the enhancement of the SdHO amphtude at 
the maximum bandwidth conditions. The interpretation 
will be confirmed in the next subsection by comparing 
the calculated traces with our experimental SdHO. 

The numerical calculation of p^lf by Shi et al7 shows 
only negligibly small share of the SdHO. This can be in- 
terpreted in terms of the factor Ji(2ttvb)/'2ttvb in eq. 
(20), keeping only the k — 1 term in the summation. 
As shown in Fig. 4 (b), the function Ji{27tvb)/2ttvb 
decrease from 0.5 with the increase of \vb\, becomes 
zero at It;^! ~ 5/8, and oscillates thereafter with ever 
decreasing amplitude of less than 0.07, with zeros at 
\vb\ ^ l/S + n/2 [n = 2,3,4,...). The oscillation of 
with \/B works just to slightly counteract the effect of 
the CO while \vb\ remains much smaller than 5/8, which 
is the case for our experiment for higher magnetic fields. 
When the amplitude of the periodic modulation is large, 
as is the case in Shi et ai, the factors cos^{qRc — 7r/4) 
and I Ji(27rws)/27ruB| conspire to alternatingly become 
small and keep the diffusion contribution to the SdHO 
small over the whole range of magnetic fields. 

3.3.2 Comparison with experimental data 

In this subsection, we compare our experimentally ob- 
tained SdHO with calculated resistivities, in an attempt 
to gain more quantitative understanding of the magnetic- 
field dependence of the SdHO amplitude. Main focus is 
on the behavior at higher magnetic-field side that re- 
mains unexplained in §3.2. 

Figures 5 and 7 show the experimentally obtained 
SdHO for Sample A and Sample B, respectively, plotted 
against 1/B. The oscillatory parts, Apsdii/po, are ob- 
tained by applying a Fourier high-pass filter to the plot 
of Pxx/po vs. 1/B. Simultaneously plotted \Vb\ serves as 
a guide to review the transition with the increase of B 
from suppression to enhancement of the SdHO ampli- 
tudes at the maxima of the band width. 

In Figs. 6 and 8, we plot calculated coUisional (hop- 
ping) and diffusion (band) contributions and the sum 



V 
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Fig. 6. (Color online) Calculated SdHO for Sample A. The col- 
lisional (hopping) contribution [eq. (21)], the diffusion (band) 
contribution [eq. (22)], and the addition of the two contributions 
are plotted (left axis), with the former two traces negatively off- 
set for clarity. The half width of the Landau bands is also plotted 
(right axis). 

of the two contributions for Sample A and Sample B, 
respectively, using corresponding sample parameters. In 
pursuit of better agreement with the experiment, we used 
slightly modified version of the argument presented in the 
previous subsection. 

As mentioned earlier, the coUisional contribution given 
by eq. (19) possesses an additional B-dependent factor 
(1 -j- 2-Kk"f) compared with the DOS in eq. (9). In prin- 
ciple, eq. (19) should also describe the SdHO of the 
plain 2DEG by placing Vq =0 {vb = 0). This, how- 
ever, is at variance with the firmly established relation 
A/9sdH/[po^(T/Tc)] (X ADi/Do^^ owing to the extra fac- 
tor, which may possibly be resulting from the approxima- 
tion used in the course of deducing eq. (12). We, there- 
fore, discard eq. (19) and assume that the relation 

Ap^ ^ Pf^ - Po ^ ^ Ai^ ^21) 
Po Po Do 

confirmed for plain 2DEGs also represents the oscilla- 
tory part of the coUisional contribution of the ULSLs. 
Note, however, that this choice does not have a drastic 
effect, since (1 -I- 2nk^) ~ 1 for large enough magnetic 
fields. In eq. (21), we assumed A{T/Tc) ~1 appropriate 
for low temperatures. The constant C has been shown 
to be equal to 2 for ideally uniform 2DEGs but deviate 
from this ideal value by a small (typically a few percent) 
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Fig. 7. (Color online) Experimentally obtained SdHO versus 1/B 
for sample B (left axis) and the half width of Landau bands (right 
axis). 
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(Color online) Similar to Fig. 6 for Sample B. 



inhomogeneity in the electron density. We selected C = 
0.92 and 2 for Sample A and Sample B, respectively, the 
values that quantitatively describe the SdHO of the ad- 
jacent plain 2DEGs and also the SdHO of the ULSLs for 
the low magnetic-field range. 

In a previous publication,^^ we have pointed out that 
a factor ^(vr/cjcTw) describing the damping of the CO 
due to scattering needs to be incorporated for quanti- 
tative account of the experimental CO traces, with Tw 



a characteristic scattering time usually identifiable with 
TQ. This will also affect the second term in eq. (20) that 
includes the CO as a multiplying factor, resulting in the 
diffusion contribution to the SdHO (excluding the first 
term corresponding to the ordinary CO) as 



,dif 
yv 
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Po 



(^cTw / akpEpTo 



{qRc 



E' 



X > ^ cos 



27rfc 



Jl{2TTkvB) ^_2^k-y 

2TrkvB 



(22) 



Thermal damping is neglected here again. The exponen- 
tial factor exp(— 277/07) still works in favor of smaller 
k. However, due to the extra i?-linear dependence men- 
tioned earlier, the more weight is on the higher magnetic 
field side for the diffusion contribution, where the ex- 
ponential factors still remain rather large. Therefore, in 
Figs. 6 and 8, we preserved up to fc = 5 in the summation 
of eq. (22). 

By comparing Figs. 5, 7 and Figs. 6, 8, it can be seen 
that the addition of the two types of contributions qual- 
itatively reproduce the experimental SdHO, notably the 
transition from the suppression to the enhancement of 
SdHO at maximum bandwidth conditions. The transi- 
tion is attributable to the rapid growth of the diffusion 
contribution with increasing B. The diffusion contribu- 
tion plays more important role in Sample B than in Sam- 
ple A. This can be ascribed to the higher mobility fi 
for Sample B; as can be seen in eq. (22), the diffusion 
contribution is proportional to t/Tq cx /i^. The effect 
of smaller modulation amplitude Vq is two fold: on one 
hand, the smaller Vq makes the counteracting effect of 
Ji{2'Kv b) / {'i.'Kv b) mentioned at the end of the preceding 
subsection less effective; on the other hand, the smaller 
Vq is disadvantageous because of the factor Vb^ in eq. 
(22). These two effects somewhat compensate each other 
to make the difference in Vb between the two samples less 
important. 

The degree of agreement between the calculated and 
the experimental SdHO shown here should be assessed 
with care, since we had to employ slightly larger values 
of Vq in eq. (21) than in eq. (22) in order to achieve 
good agreement with experimental traces, as mentioned 
in §3.2. Although the inconsistency should be resolved 
in the future studies, it does not affect the qualitative 
argument presented here. 

4. Conclusions 

We have shown that the experimentally observed 
SdHO for a ULSL is basically reproduced by the addi- 
tion of the coUisional contribution Apl^"' [eq. (21) with 
eq. (10)] and the diffusion contribution Ap|Jl^ [eq. (22)]. 
The amplitude of the oscillation alternates between the 
two contributions: Ap™' is suppressed while ApJ^|^ is en- 
hanced at the maximum bandwidth conditions [eq. (4)]. 
Owing to an extra linear- _B factor for ApJ^i'^ in addition 
to the common exponential damping factor, Ap™' dom- 
inates the SdHO at low magnetic fields but ApJ^J^ out- 
weighs Ap™' at higher magnetic fields. This accounts for 
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the experimentally observed transition from suppression 
to enhancement of the SdHO at the maximum bandwidth 
conditions. The term Ji{2ttvb) / {"^t^vb) in eq. (22) quali- 
tatively explains why diffusion contribution to the SdHO 
was not observed in a previous experiment.^ 
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Appendix: Derivation of eq. (16) 

In this appendix, we describe the derivation of eq. (16) 
from eq. (14). Performing the summation first, eq. (14) 
becomes 



Performing the integral in eq. (A-1), we finally achieve. 
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where we introduced the notation a = ef — I/^ — vb cos t 
for brevity. Noting that 2A^ + 1 = 2a + 1 + 2(Af - a), the 
summation in the last line of eq. (A-2) can be performed 
by appealing to the approximation J2n=o ~^ J2n=-oo 
as before, and by applying the Poisson sum formulae, 
resulting in 
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Replacing this result into eq. (A-2), we obtain 
S{EF,t)= ^ 
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Owing to the smallncss of To /Ep and Vb / E-p , the first 
term is by far the dominant term, ending up with eq. 
(16). 
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approximation exerts negligible effect on D{E) for F E. 
Here, the Ey is approximated by a constant value at B = 
0. Strictly speaking. Bp oscillates with B to fulfill rie = 
f{E)D{E)dE, with f{E) the Fermi-Dirac distribution 
function. However the effect is small at low magnetic fields 
and is negligible in the context of the present study. 
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